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A phenomenological explanation of the magnetoelectric behavior of CosTeOe is developed. We 
explain the second harmonic generation data and the magnetic field induced spontaneous electric 
polarization in the magnetically ordered phase below 20K. 
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I. INTRODUCTION 

Recently there has been an explosion in the number of 
compounds which exhibit nontrivial magnetoelectric be- 
havior at low temperatures. Co 3 Te0 6 (CTO) is an 
interesting such system whose properties have recently 
been studied. [6j Although the magnetic structure is as 
yet not clarified, it seems useful to construct a mean- 
field scenario which can explain the major experimen- 
tal results. The measurments of Hudl et al.[6| of M/H, 
d(M/H)/dT, and C/T versus T, where M is the magne- 
tization, H the magnetic field, and C the specific heat, 
indicate that there are at least two magnetic phase transi- 
tions at temperatures below about 30K, one at T\ « 26K 
and another at « 18. 5K, but the details of the mag- 
netic tructure are not known, other than that the system 
is not ferromagnetic. According to Ref. 0, the mag- 
netic structure is described by several incommensurate 
wave vectors. Single crystal netron diffraction measure- 
ments reveal that the incooomensurate wave vector (s) are 
in the a-b plane and not along c.Q We propose the ex- 
istence of magnetic order at zero wave vector, consistent 
with the results of Li et ai, although, this would have 
to involve an antiferromagnetic arrangement of moments 
within the unit cell to give the observed zero net mo- 
ment. In addition, our analysis suggests the appearance 
of an additional magnetic phase transition. In the ab- 
sence of magnetic order the crystal symmetry is [Hill that 
of space group C2/c (#15 in Ref. [To]). We will take the 
generators of this space group to be the glide operation 
TOft = (x, —y, 2+1/2), a two-fold screw rotation about the 
crystal b axis, 2;, = (— x, y + 1/2, — z+1/2), and the three 
translations, (a: + 1/2, y + 1/2, z), (x - 1/2, y + 1/2, z), 
and (x, y, z + 1), where x, y, and z are in units of lattice 
constants. These are equivalent to those of Ref. [ll|. 



II. EXPERIMENTAL DATA 

The data of Hudl et aL@ consist of several types. 
As mentioned above, the measurements of magnetiza- 
tion and specific heat indicate phase transitions at at 
least two temperatures, T± and T%, but the nature of 
magnetic ordering could not be determined from their 
data. The lower-temperature transition may be a dis- 
continuous one. Of primary interest to us is their mea- 



surement of the intensity of second harmonic generation 
(SHG), whose cross section is proportional to the third 
order electric susceptibility x Q ^ 7 , where a, (3 and 7 label 
components (or in the present case label crystallographic 
directions). Their experimental geometries are chosen 
such that the SHG cross section is proportional to Xaaa ■ 
In a system having high symmetry, e. g. having inversion 
symmetry, the SHG intensity is zero for all frequencies, 
and this applies to CTO above about T2 = 18. 5K. How- 
ever, below that temperature they find that Xaaa and 
Xccc are nonzero, but Xbbb is apparently zero at all tem- 
peratures. From this they conclude that the point group 
retains only m,b symmetry. As we shall see, if, as they 
assert, the symmetry is magnetically broken, this is not 
a correct conclusion. 

Another type of data of crucial interest to us the mea- 
surement of the electric polarization, P in the a-c plane, 
as a function of temperature and magnetic field for mag- 
netic fields along the crystallographic a and c directions. 
For zero magnetic field, at temperatures below about 
18K, they find a very small, possibly zero, spontaneous 
polarization in the a and c directions which increases al- 
most proportional to the magnetic field. In fact, we find 
that their results for P c at T = 5K as a function of H a 
can be fit within experimental uncertainty (±5 in P c ) to 



P c = -0.15 + 6.93^ + 0.33^ 



(1) 



with P c in /xC/m 2 and H in Tesla. In other words, they 
found an important magnetic field-dependent contribu- 
tion to P c linear in H a with P c (H a = 0) « 0. 



III. SYMMETRY ANALYSIS 

We will carry out our analysis in terms of an expansion 
about the "vacuum" , which we take to be the phase above 
26K in which the magnetic order parameters and electric 
polarization are zero. Magnetically ordered phases are 
described by nonzero magnetic order parameters. We 
will also discuss briefly nonmagnetic structural distor- 
tions which lower the crystal symmetry from C2/c and 
which are described by appropriate order parameters. 
Although, as mentioned in Ref. [f| there may exist in- 
commensurate magnetic order described by M(q) with 
q ^ 0, incommensurate magentic order can not, by it- 
self, explain the experimental results, as we will explain 



below. 



A. Electric Polarization 

We first review the phenomenological theory of mag- 
netization induced electric polarization P. The magne- 
toelectric free energy is of the form 

Fme = ^P 2 +E Ai?(n) ' W 

n 

where \e is the dielectric susceptibility of the vaccum 
(the phase above T — 26K), which we assume to be 
isotropic for simplicity and AF^- 71 ' is the contribution lin- 
ear in P (so that it induces a nonzero value of P) and of 
order H n . For instance, to lowest order in powers of the 
magnetic order parameters, we write [5lfl2| 

AF (o) = ^a aW (q)-P«[^(q)*^(q)-^fe(q)^(q)*] 

+b akl P a M k {q = 0)M l (q = 0) 
AF^ = c af3k P a HpM k {q = 0) 
(q)P a H p H 1 

x[M fc (q)*Af,(q) - M fc (q)Mj(q)*] 
+e a p jk iP a H p H y M k (q = 0)M l (q = 0) , (3) 

where we invoke the Einstein convention which implies 
summation over repeated subscripts, Greek subscripts la- 
bel crystallographic directions, and Roman letters label 
irreducible representations (irreps), which in the present 
case are one dimensional. The magnetic order parameter 
Mfc(q) can be thought of as the amplitude of the mag- 
netic normal mode associated with irrep . 5] These nor- 
mal modes are the linear combinations of magnetic mo- 
ments within the unit which bring the quadratic terms 
in the Landau expansion into diagonal form. We will 
discuss the symmetry of the M k s in a moment. Here 
the Fourier transforms are defined so that for q ^ 0, 
IMfc(q) = Mfe(q)*, where I = m b 2 b is spatial inversion. 

-Fme must be invariant under all the symmetries of the 
"vacuum." These symmetries include time reversal sym- 
metry, transnational symmetry (which leads to wave 
vector conservation), and the crystallographic symme- 
tries m b and 2 b (which together imply invariance under 
spatial inversion I). We will consider the crystallogrphic 
symmetries in a moment. Time reversal symmetry re- 
quires that the total number of powers of H and M(q) 
must be even. The condition that Fme be real valued 
implies that a(q) and d(q) be pure imaginary. The form 
of AF^ 1 ' is such that wave vector conservation implies 
that the magnetic order for this mechanism must occur 
at zero wave vector, and, as previously noted, it must 
be antiferromagnetic to be consistent with the observed 
zero net magnetic moment of the system. (In fact CTO 
has a large enough paramagnetic unit cell that antiferro- 
magnetic order can develop without increasing the size 



of the unit cell, as occurs in LaTiOaflU and C^Osflij.) 
Such an antiferromagnetic moment would be consistent 
with the magnetic measurements of Hudl et al. [6\ 

When Fme is minimized with respect to P to obtain 
its equilibrium value, one sees that AF^ gives rise to a 
contribution to P which is of order H n . In many mul- 
tiferroics, such as Ni 3 V 2 8 @|3 (NVO) and TbMnO 3 
(TMO), AF^ is a crucial term which gives rise to a 
spontaneous polarization at H = 0. Many other cases 
are similarly analyzed in Ref. 0. In these cases, the mag- 
netic order is incommensurate, so that the polarization 
(a zero wave vector property) can not be linear in the 
magnetic order parameter. Since in CTO P oc H, we 
consider AF^ from which we get 

P a = XEY Ca l 3kH l 3Mk ■ ( 4 ) 

We now show how the crystallographic symmetries 
constrain the coefficient tensor c a p k . In particular, we 
will show that these symmetries fix the symmetry of M k ■ 
For this purpose, note that AF^ has to be invariant 
under these symmetries. In this analysis, we will confine 
P and H to be perpendicular to the crystallographic b 
direction, as they were in the experiments of Ref. |6|. In 
that case, we only consider terms in Ai" 1 ' with a and 
ft labeling the crystallographic a and c directions, and k 
labels the possible magnetic irreps at zero wave vector. 
Remembering that H is a pseudovector, we note that 

m b [P a H ] = -P a Hp , 2 b [P a H p ] = P a H p .(5) 

Accordingly, for AF^ to be an invariant we require that 

m b M k = -M k , 2 b M k = M k . (6) 

To implement Eq. ([6]) we need to characterize the sym- 
metry of the magnetic ordering, which we have inferred 
occurs at zero wave vector. For phase transitions the cat- 
alog of broken symmetry phases that can result from a 
phase transition in any of the 230 crystallographic space 
groups can be obtained using the suite of computer pro- 
grams ISODISTORT which is accessible on the web.(l5j 
As applied to CTO one predicts that only four magnetic 
irreps can result from a single phase transition at zero 
wave vector. This formulation specifically does not allow 
for a multicritical point at which there is a simultaneous 
breaking of two distinct symmetries. For CTO there is 
no experimental indication that the magnetic phase tran- 
sitions arise from such a multicritical point. fill] There- 
fore we assume the validity of the four possible magnetic 
phases of Table [fl which ISODISTORT lists for space 
group C2/c. Looking at Table U we see that to be con- 
sistent with Eq. |6]), the magnetic order parameter can 
only be that of irrep T2- 

B. Second Harmonic Generation 

We now turn to the analysis of the SHG cross sec- 
tion at H = 0. To develop a nonzero SHG cross sec- 
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TABLE I: Symmetry of the magnetic irreps T„ at zero wave 
vector for CTO. Here X(O) is the eigenvalue of the operator O: 
OM k = \{<D)M k , where M k = M(T k ) is the order parameter 
associated with the fcth irrep. Also £ is the identity and T is 
spatial inversion. In the last line, we give the direction of the 
ferromagnetic moment if it is allowed to be nonzero. 
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tion a quantity like dxaaa/dMp must be nonzero in the 
vacuum (magnetically disordered phase), so that when 
we turn on the magnetic order parameter Mp (in the 
magnetically ordered phase) the SHG cross section be- 
comes nonzero. To study this quantity it is useful to 
note that it has the symmetry of d[p a p a p a }/dMk, where 
p a is the a-component of the dipole moment operator 
and Mk is a magnetic order parameter. One sees that 
this quantity is zero because M& is odd under time re- 
versal, and the dipole moment operator is even under 
time reversal. [TtJ Therefore, the phenomenological expla- 
nation for a nonzero SHG cross section must come from 
X a = d 2 {p a p a p a ]/[dM k (q)dMj*(q)} being nonzero in the 
disordered phase. This quantity has the same symmetry 
as X a ee p\M where M = M fe (q)Mf (q) or M = M k M h 
The fact that the SHG is proportional to the product of 
two different order parameters, each of which, as we shall 
see, describes a one dimensional irrep, has been noted 
before. 18]. Here, from the polarization data, we know 
of the existence of at least one irrep at zero wave vector 
and according to Ref. magnetic ordering occurs with at 
least one irrep at nonzero wave vector. To have a nonzero 
SHG cross section we need a second irrep, either at zero 
wave vector or at the same nonzero wave vector. In either 
case the appearance of a second irrep requires an as yet 
unobserved phase transition, which may be unobtrusive 
enough that it was not seen by Hudl et at. We consider 
these two scenarios in turn. 

The condition for a nonzero SHG cross section is iden- 
tical to that for a nonzero electric polarization because 
the symmetry properties of the dipole moment operator 
and the electric polarization are the same. Thus, if Xaaa 
and Xccc are nonzero, then P a and P c are expected to 
be nonzero. Furthermore, no matter which scenario is 
adopted, there is a possible problem in that although ex- 
periments show that for H — 0, Xaaa and Xccc are nonzero 
and Xbbb — 0, the expected field independent contribu- 
tions to P a and P c are very small. The explanation for 
this may be that the SHG is anomalously large when the 
polarization is due to modification of electronic orbits (as 
contrasted to being due to ionic displacements) . plj 

In the first scenario, we assume that the nonzero 
SHG cross section is induced by magnetic order at zero 



wave vector and study the symmetry properties of X a . 
Since p 2 a transforms like unity, it suffices to study X a = 
p a MkMi, to indicate whether Xaaa is or is not zero. 
Since Xbbb = 0, we require that pbM k Mi be odd under ei- 
ther TOb or 2b. This implies that M k Mi either be even un- 
der nib or odd under 2b. Using Table HT1 we see that this 
criterion excludes either M\M 2 or M3M4 being nonzero. 
Similarly if Xaaa and x.ccc are nonzero, we require that 
both p a MkMi and p c M k Mi be even under both m& and 
2b. This implies that M k M\ be even under m& and odd 
under 2b. These requirements indicate that either M1M4 
or M2M3 be nonzero. Since we have previously invoked 
the existence of irrep M2 to explain the electric polar- 
ization, we opt for M2M3 being nonzero. The fact that 
the magnetic moment perpendicular to b (coming from 
irrep M3) is zero (or very small) would have to be a result 
specific to the details of the interactions. 

In the second scenario one would have to posit an ad- 
ditional phase transition involving a second incommen- 
surate magnetic irrep to give rise to a nonzero SHG cross 
section. In principle, one would have an accompany- 
ing field independent polarization coming from Af<°) 
whose absence in experiment would have to be explained 
as above in terms an unusually large SHG cross section. 
To illustrate this mechanism consider the hypothetical 
case when the incommensurate magnetic ordering occurs 
at q = qob. In this case one finds that there are two mag- 
netic irreps, one of which, call it Mi(q), is even under 2& 
and the other, call it M2(q), is odd under 2b. Then 
one sees that X = p [Mi(q)*M 2 (q) - Mi(q)M 2 (q)*] 
and Y = p c [Mi(q)*M 2 (q) - Mi(q)M 2 (q)*] are both in- 
variant under 2b (and under Z), so that Xaaa X and 
Xccc Y are allowed to be nonzero, whereas Xbbb remains 
zero. In a common scenario Q one irrep would give rise 
to nonzero magnetic moments along the b axis, and the 
other would give rise to nonzero magnetic moments along 
the c axis. These irreps would be out of phase (so that 
Afi(q)*M 2 (q) - MjJa)M 2 (q)* is nonzero) giving rise to 
a magnetic spiral. [19[ 

TABLE II: As Table U Symmetry of the product of two zero 
wave vector magnetic irreps F n for CTO. 
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C. Discussion 

To summarize our conclusions: we require the exis- 
tence of zero wavevector magnetism according to irrep 
M 2 to explain the magnetic field induced electric polar- 
ization. In one scenario we explain the SHG cross section 
as being proportional to M2M3. Since we prefer not to 
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assume a multicritical point, the latter result would im- 
ply that there are actually two phase transitions. At the 
higher-temperature transition (at T = 18.5K) a magnetic 
field induced spontaneous electric polarization appears 
and at the lower-temperature transition (at some tem- 
perature close to but below 18. 5K) the SHG cross section 
becomes nonzero. Here a very small magnetic field inde- 
pendent polarization should also appear. In principle, 
one would hope to show the temperature dependence of 
the SHG cross section to be proportional to the product 
of these two order parameters whose temperature depen- 
dence was independently established by neutron diffrac- 
tion. This type of experimental program was carried out 
for the electric polarization of NVO (see Fig. 6 of Ref. 
I20I ). Note also a magnetically induced SHG cross section 
implies that the symmetry involves time reversal. The 
magnetic phase with irrep Mi is odd under mi,, as indi- 
cated in Table |U In contrast, if we were dealing with a 
nonmagnetic structural phase transition, as the analysis 
of Hudl et al. tacitly assumes, then the low-temperature 
phase would be even under nib, as they state. However, 
note that the presence of magnetic irreps M2 and M3 
breaks the mirror symmetry of mb, but the symmetry 
of mi, plus time reversal is maintained. This is consis- 
tent with the results of Tables 7 and 4 of Ref. [H (The 
misidentification of Ref. @ is not completely harmless. 
If one assumes that mj symmetry is unbroken, then, as 
they find, it is impossible to use AF^ to explain why 
dPa/dHp is nonzero for a,/3 = a, c.) 

The second scenario has similar ramifications except 
that it involves magnetic ordering at some incommen- 
surate wave vector. This scenario would also require a 



second phase transition at which a second incommensu- 
rate order parameter would appear. [23| In principle, such 
a transition could involve a slightly different wave vector 
than that already present. But, as argued in Ref. 0, quar- 
tic terms in the Landau free energy would favor locking 
these two nearby wave vectors to the same value. 

We have implicitly assumed that the experimental re- 
sults are induced by magnetic ordering. One might ques- 
tion whether the results of Ref. could be explained by 
simply invoking one or more phase transitions driven by 
structural distortions. Since magnetic ordering appears 
at these transitions the question is which order parame- 
ter is the primary one whose presence induces the appear- 
ance of the other one. If Q is a structural order parameter 
(like the tilting angle of a cage of oxygen ions), then one 
can invoke an interaction of the type V ~ M(Tk)M(Ti)Q 
to explain the appearance of a nonzero value of Q at the 
transition. Via this coupling the appearance of one or 
more magnetic order parameters (which are the primary 
order parameters) would induce a structural distortion 
(because Q appears linearly). The converse case, where 
the magnetic order parameter appears linearly and the 
primary order parameter Q appears quadratically (or lin- 
early, for that matter) is not allowed by time reversal 
symmetry. But if the magnetic order parameters are 
the primary ones, then the theoretical approach of the 
present paper is essentially unchanged by the appearance 
of secondary structural order parameters. 
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